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Abstract 



We discuss a minimal Supersymmetric SO(10) model where B-L symmetry is 
broken by a 126 dimensional Higgs multiplet which also contributes to fermion 
masses in conjunction with a 10 dimensional superfield. This minimal Higgs 
choice provides a partial unification of neutrino flavor structure with that of 
quarks and has been shown to predict all three neutrino mixing angles and 
the solar mass splitting in agreement with observations, provided one uses 
the type II seesaw formula for neutrino masses. In this paper we generalize 
this analysis to include arbitrary CP phases in couplings and vevs. We find 
that (i) the predictions for neutrino mixings are similar with U e3 ~ 0.18 as 
before and other parameters in a somewhat bigger range and (ii) that to 
first order in the quark mixing parameter A (the Cabibbo angle), the leptonic 
mixing matrix is CP conserving. We also find that in the absence of any higher 
dimensional contributions to fermion masses, the CKM phase is different from 
that of the standard model implying that there must be new contributions 
to quark CP violation from the supersymmetry breaking sector. Inclusion of 
higher dimensional terms however allows the standard model CKM phase to 
be maintained. 



I. INTRODUCTION 

The observation of solar and atmospheric neutrino deficits in various experiments such as 
Chlorine, Super-Kamiokande, Gallex, SAGE and SNO together with recent results from the 
K2K and Kamland [1] experiments that involve terrestrial neutrinos have now conclusively 
established that neutrinos have mass and they mix among themselves. In conjunction with 
negative results from CHOOZ and PALO- VERDE reactor experiments, one now seems to 
have a clear idea about the mixing pattern among the three generations of neutrinos. Of 
the three angles needed to characterize these mixings i.e. 12 , #23 and 6*i 3 , the first two are 
large and the third is small [2]. 

As far as the mass pattern is concerned, there are three distinct possibilities: (i) "normal 
hierarchy" where mi <C m 2 <C m 3 ~ \J Am\; (ii) inverted one where m\ ~ m 2 — J Am\ ^> 
m 3 ; (iii) quasi-degenerate, where mi ~ m 2 — m 3 . 
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The normal hierarchy is similar to what is observed in the quark sector although the 
degree of hierarchy for neutrinos is much less pronounced. Understanding this difference is 
somewhat of a challenge although grand unified theories generally imply similar hierarchies 
for quarks as well as neutrinos and therefore provide a partial explanation of this challenge. 

The observed mixing pattern for leptons however, is totally different from what is ob- 
served for quarks, posing a much more severe theoretical challenge, in particular for grand 
unified theories that unify quarks and leptons. Both these questions have been the subject 
of many papers [3] that propose different approaches to solve these problems. Our goal in 
this paper is to address them in the context of a minimal SO (10) model, which provides an 
interesting way to resolve both the mass and mixing problems without making any extra 
assumptions, other than what is needed to derive the supersymmetric standard model from 
SO(10). 

The first question one may ask is: why SO(10) ? The answer is that we seek a theory 
of neutrino mixings that must be part of a larger framework that explains other puzzles of 
the standard model such as the origin of matter, the gauge hierarchy problem, dark matter 
etc. As is explained below, supersymmetric SO(10) has the right features to satisfy all the 
above requirements. 

Supersymmetric extension of the standard model (MSSM) seems to be a very natural 
framework to solve both the gauge hierarchy as well as the dark matter problems. The 
nonrenormalization theorem for the superpotential solves an important aspect of the gauge 
gauge hierarchy problem i.e. the radiative corrections do not destabilize the weak scale. As 
far as the dark matter puzzle goes, if one adds the extra assumption that R-parity defined 
as R = (— 1) 3 (- B_i ) +2,S ' is also a good symmetry of MSSM, then the lightest supersymmetric 
particle, the neutralino is stable and can play the role of dark matter. In fact currently all 
dark matter search experiments are geared to finding the neutralino. An added bonus of 
MSSM is that it also provides a nice way to understand the origin of electroweak symmetry 
breaking. 

Turning to neutrinos, the seesaw mechanism [4] which is the simplest way to understand 
small neutrino masses, requires a right handed neutrino with a large Majorana mass Mr. 
Present data require that M R has to be very high and yet much smaller than the Planck mass. 
This would suggest that perhaps there is a new symmetry whose breaking is responsible 
for M R ; since symmetry breaking scales do not receive destabilizing radiative corrections 
in supersymmetric theories, one has a plausible explanation of why M R <C Mp£. One 
such symmetry is local B-L symmetry [5], whose breaking scale vbl could be responsible 
for M R (i.e. M R = fv B L) since the righthanded neutrino has nonzero B-L charge and is 
massless in the symmetry limit. The formula for neutrino masses in the seesaw model is 
given by m v ~ — j^', therefore their smallness is now very easily understood as a result 
of mo <C vbl 1 - The smallest grand unification group that incorporates both the features 
required for the seesaw mechanism i.e. a righthanded neutrino and the local B-L symmetry 
is SO(10). 



The symmetry could also be a global symmetry without conflicting with what is known at low 
energies [6]; but we do not discuss it here. 
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Secondly, the seesaw mechanism also goes well with the idea of grand unification. It 
is well known that if MSSM remains the only theory till a very high scale, then the three 
gauge couplings measured at LEP and SLC at the weak scale can unify to a single coupling 
around M v ~2x 10 16 GeV. Curiously enough, naive estimates for the seesaw scale needed 
to understand the atmospheric neutrino oscillation within a quark-lepton unified picture 
implies that the seesaw scale must be around vbl > 10 15 GeV. Close proximity of vbl and 
Mjj suggests that one should seek an understanding of the neutrino masses and mixings 
within a supersymmetric grand unified theory. While this already makes the case for SUSY 
SO (10) [7] as a candidate theory for neutrino masses quite a strong one, there is an additional 
feature that make the case even stronger. 

This has to do with understanding a stable dark matter naturally in MSSM. From the 
definition of R-parity given above, it is clear that in the limit of exact SO(10) symmetry, 
R-parity is conserved [8]. However, since ultimately, SO (10) must break down to MSSM, 
one has to investigate whether in this process R-parity remains intact or breaks down. This 
is intimately connected with how the right handed neutrino masses arise (or how local B-L 
subgroup of SO(10) is broken). It turns out that the B-L symmetry needs to be broken by 
126 Higgs representation [9-11] rather than 16 Higgses [12] as explained below. 

It was pointed out in [9] that SO(10) with one 126 and one 10 is very predictive in 
the neutrino sector without any extra assumptions, while at the same time correcting a 
bad mass relation among charged fermions predicted by the minimal SU(5) model, i.e. 
m^Mu) = m s (Mjj). There is a partial unification of flavor structure between the neutrino 
and quark sectors and since all Yukawa coupling parameters of the model are determined 
by charged fermion masses and mixings, there are no free parameters (besides an overall 
scale) in the neutrino mass matrix. It is therefore a priori not at all clear that the neutrino 
parameters predicted by the model would agree with observations. In fact the initial analysis 
of neutrino mixings in this model that used the type I seesaw formula and did not include CP 
phases did predict neutrino mixings that are now in disagreement with data. In subsequent 
papers [13,10,14-16], this idea was analyzed ( in some cases by including more than one 10 
Higgses) to see how close one can come to observations. The conclusion now appears to be 
that one needs CP violating phases to achieve this goal [16] 2 . 

Another approach is to use the type II seesaw formula for neutrino masses [18] i.e. 

M v ~ f-^h -J^—. (1) 

In models which have asymptotic parity symmetry such as left-right or SO(10) models, it is 
the type II seesaw that is more generic. 

A very interesting point about this approach, noted recently [19], is that use of the type 
II seesaw formula for the two generation subsector of and u T , and dominance of the first 
term leads to a very natural understanding of maximal atmospheric mixing angle due to 
b — r mass convergence at high scale. 

Whether the above idea does indeed lead to a realistic picture for all three neutrino 
generations was left unanswered in ref. [19]. To answer this question, a complete three 



2 For a different class of SO(10) models with 126 Higgs fields see [17] 
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generation analysis of this model was carried out in [20] where it was pointed out that the 
same b — r convergence condition that led to large atmospheric mixing angle also leads to 
a large solar angle and also a small 9 13 . Furthermore, it also resolves the mass puzzle for 
neutrinos (i.e. a milder hierarchy for neutrinos than that for quarks) since it predicts that 



y Am| / Am^ ~ A ~ 0.22, where A is the small quark mixing parameter in the Wolfenstein 

parametrization of the CKM matrix. As there are no free Yukawa coupling parameters in 
the neutrino sector, it is quite amazing that all the neutrino parameters can come out in 
the right range. 

In the analysis of ref . [20] , CP violating phases were set to zero. It was implicitly assumed 
that all known CP violating processes in this model would arise from the supersymmetry 
breaking sector which would then make the model completely realistic even though it does 
not have the conventional CKM CP violation. CP violation is however a fundamental 
problem in particle physics and its origin at the moment is unclear. It is therefore of interest 
to see (i) whether the minimal SUSY SO(10) can remain predictive for neutrinos even after 
the parameters in the model are allowed to become complex, thereby ushering in CP violation 
into the quark sector in a direct way and (ii) if any other useful information on the nature 
of CP violation for both quarks and leptons can be gained in this model. 

To study CP violation in this model, we generalize our earlier analysis to make all 
Yukawa coupling parameters as well as the vacuum expectation values (vev) complex. One 
might suspect that since this will bring in several new parameters, the model may lose its 
predictivity. We find that this is not so. We have seven phases, one of which can play the 
role of the CKM phase. Due to the presence of a sum rule involving the charged lepton and 
quark masses, despite the presence of extra phases, the model still remains predictive and 
pretty much leads to the same predictions with minor changes as in the case without phase. 

Assuming that b-tau mass convergence leads to maximal neutrino mixing, constrains 
three of the phases to be equal and matching the electron mass fixes two others. The 
remaining arbitrary phase is associated with the up quark, whose tiny mass keeps this phase 
well hidden from this discussion. The CKM phase however turns out to be outside the one a 
region of the present central value in the standard model fits [21]. That means that one will 
need some contributions to the observed CP violation from the SUSY breaking sector. We 
then observe that if we include the higher dimensional contributions to the fermion masses 
which were ignored before, only for the first generation (which can be done naturally using 
an R-symmetry), this introduces only one new parameter which relaxes the electron mass 
constraint but does not affect the neutrino sector. We can maintain a CKM phase equal to 
that given by the standard model fit (of about 60°) . 

An important outcome of this analysis is that in both cases, to first order in Cabibbo 
angle A, the leptonic mixing matrix is CP conserving, which can therefore be used to test 
the model. 

In our opinion, these observations have lifted the minimal SO(10) with 126 to a realistic 
grand unification model for all forces and matter and ought therefore be considered as a 
serious candidate for physics beyond the standard model up to the scale of 10 16 GeV. Just 
like the minimal SUSY SU(5) grand unified theory could be tested in proton decay searches, 
this minimal version of SO(10) can be tested by future neutrino experiments. Important 
experiments for this purpose are the planned long base line experiments which will provide 
a high precision measurement of the mixing parameter 9 13 (also called U e3 ) to the level of 
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0.1 or less. 

This paper is organized as follows: we start in sec. 2 with a few introductory remarks 
about minimal SO(10) with 126 Higgses; we then sketch a derivation of the type II seesaw 
formula in sec. 3 and show how b — r convergence leads to an understanding of large neutrino 
mixings; in sec 4, we discuss the minimal SO(10) model with all couplings and vevs real and 
discuss the prediction for neutrinos. Section 4 gives more details of our earlier paper [20] 
and gives an analytic explanation of the constraints on the model parameters. In sec. 5, we 
discuss the effects of including the most general form of CP violating phases and present our 
predictions for neutrino masses and mixings without the higher dimensional operators; in 
this section, we also comment on the implication of including higher dimensional operators. 
In sec. 7, we summarize our results and present our conclusions. 

II. A MINIMAL SO(10) MODEL 

In any SO(10) model, one needs several multiplets to break the symmetry down to 
SU(3) C x U{l) em . Usually, to break the SO(10) group down to the Pati-Salam group 
SU(2) L x SU(2) R x SU{A) C or to the left-right group SU(2) L x SU(2) R x U(1) B -l x SU(3) c , 
one needs either 54 or 54©45 Higgs multiplets [10]. Then to break SU(2) R x U(1) B -l 
down to U(1)y, one can employ either 16 + 16 or 126 + 126 pair since both these models 
have standard model singlet and B-L^O fields in them. The reason for having the complex 
conjugate is to maintain supersymmetry down to the electroweak scale. Finally to break the 
standard model group, 10 is used. Thus one needs always at least five Higgs multiplets in 
most constructions of SO (10) model. One could replace 54+45 pair by a 210 representa- 
tion [22], reducing the number of multiplets required to four. The neutrino results that we 
discuss are not affected by the choice of Higgses that affect the breaking in the first stage 
but are crucially dependent on how one implements the subsequent ones. 

In our model, we will use 126 + 126 to break B-L symmetry for the following reasons. 

A. R-parity and 16 vs. 126 

As discussed in the introduction, an important argument in favor of MSSM being the 
TeV scale theory is the possibility that the lightest SUSY partner can play the role of dark 
matter. In fact a lot of resources are being devoted to discover the supersymmetric dark 
matter particle. For MSSM to provide such a dark matter particle, it is important that it has 
R-parity conservation. The MSSM by itself does not have R-parity and ad hoc symmetries 
are stuck into the MSSM to guarantee the existence of stable dark matter. SO(10) provides 
an interesting way to guarantee automatic R-parity conservation without invoking any ad 
hoc symmetry as we see below. 

The crucial question is whether B-L subgroup is broken by a (i) 16-dimensional Higgs 
field [12] or (ii) 126-dimensional ones [9-11]. 

In case (i), B-L symmetry is broken by a Higgs field (the z/R-like Higgs field in 16) which 
has B — L = 1. If one looks at higher dimensional contributions to the superpotential of the 
form 16^16#, in terms of the MSSM superfields they have the form QLd c v c Hl u c d c d c v c H etc 
where Q, L, u c , d c , e c are the matter superfields of the MSSM. The field is the Higgs field 
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in 16 that breaks B-L via < v c H >= vbl- After symmetry breaking, these nonrenormalizable 
couplings will induce R-parity breaking terms of MSSM such as QLd c and u c d c d c etc with 
a slightly suppressed coupling i.e. vbl/M P £ ~ 1CT 2 — 1CT 3 . This suppression is not enough 
to let the neutralino play the role of cold dark matter - not to mention the fact that it leads 
to extremely rapid proton decay. 

It is sometimes argued that the final theory from which this effective SO(10) model 
emerges may have additional local U(l)x symmetries that will prevent these dangerous 
higher dimensional terms. To see what this implies, one can imagine a Higgs field X which 
is charged under U(l) x and has vev of order of the GUT scale. The leading order operator 
that will keep the theory safe from proton decay problem has to involve an operator of the 

form u c d c d c (jj^-) and similarly for other operators. This means that the U(l) x charge of 
X and the operator u c d c d c must be arranged in a very specific way, which leads to another 
kind of naturalness problem. 

On the other hand if B-L symmetry of SO(10) is broken by a 126 Higgs field, as in case 
(ii), the B-L^O and standard model singlet field that breaks B-L has B-L=2. As a result 
after this Higgs field acquires a vev, a Z 2 subgroup of B-L still survives and it keeps R-parity 
as a good symmetry. This has been established by a detailed analysis of the superpotential 
in Ref. [11]. This not only forbids dangerous baryon number violating terms but also allows 
for the existence of a neutralino dark matter without the need for any additional symmetry. 
We will therefore work with an SO(10) model where where the only field that breaks B-L is 
in the 126-dimensional representation. 

B. Mass sumrules in minimal SO(10) 

A second advantage of using 126 multiplet instead of 16 is that it unifies the charged 
fermion Yukawa couplings with the couplings that contribute to righthanded as well as 
lefthanded neutrino masses, as long as we do not include nonrenormalizable couplings in 
the superpotential. This can be seen as follows [9,10]: it is the set 10+126 out of which 
the MSSM Higgs doublets emerge; the later also contains the multiplets (3, 1, 10) + (1, 3, 10) 
which are responsible for not only lefthanded but also the right handed neutrino masses in 
the type II seesaw formula. We explain this below. Therefore all fermion masses in the 
model are arising from only two sets of 3 x 3 Yukawa matrices one denoting the 10 coupling 
and the other denoting 126 coupling. 

In view of the above remarks, the SO(10) model that we will work with in this paper has 
the following features: It contains three spinor 16-dim. superfields that contain the matter 
fields (denoted by ip a ); two Higgs fields, one in the 126-dim representation (denoted by A) 
that breaks the SU{2)r x U{1)b-l symmetry down to U{1)y and another in the 10-dim 
representation (H) that breaks the SU(2) L x U(l)y down to U(l) em . The original SO(10) 
model can be broken down to the left-right group SU(2) L x SU(2) R x U{\)b-l by 54 © 45 
Higgs fields denoted by S and A respectively. 

To see what this model implies for fermion masses, let us explain how the MSSM doublets 
emerge and the consequent fermion mass sumrules they lead to. As noted, the 10 and 126 
contain two (2,2,1) and (2,2,15) submultiplets (under SU(2) L x SU(2) R x SU(A) C subgroup 
of SO(10)). We denote the two pairs by and A Mjd . At the GUT scale, by some doublet- 
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triplet splitting mechanism these two pairs reduce to the MSSM Higgs pair (H u , H d ), which 
can be expressed in terms of the <fi and A as follows: 

H u = cos a u (p u + e llu sin a u A u (2) 
H d = cos a d <p d + e tJd sin a d A d 

The details of the doublet-triplet splitting mechanism that leads to the above equation are 
not relevant for what follows and we do not discuss it here. As in the case of MSSM, we will 
assume that the Higgs doublets H u>d have the vevs < H® >= v sin/3 and < H d >= vcos/3. 

In orders to discuss fermion masses in this model, we start with the SO(10) invariant 
superpotential giving the Yukawa couplings of the 16 dimensional matter spinor ipi (where 
i, j denote generations) with the Higgs fields H 10 = 10 and A = 126 3 . 

W Y = hijipiipjHw + fijipiipjA (3) 

SO (10) invariance implies that h and / are symmetric matrices. We ignore the effects coming 
from the higher dimensional operators, as we mentioned earlier. 

Below the B-L breaking (seesaw) scale, we can write the superpotential terms for the 
charged fermion Yukawa couplings as: 

W = h u QH u u c + h d QH d d c + h e LH d e c + fiH u H d (4) 

where 

h u = h cos a u + fe llu sin a u (5) 
h d = h cos a d + fe tJd sin a d 
h e = h cos a d — 3/e rw sin a d 

In general a u ^ a d and this difference is responsible for nonzero CKM mixing angles. In 
terms of the GUT scale Yukawa couplings, one can write the fermion mass matrices (defined 
as C m = ^M^r) at the seesaw scale as: 

M u 
M d 

M e 

where 

h = h* cos a u sin (3 (7) 

/ = f*e ilu sina u sin/3 
cosa d 

T\ = COt p 

coso; u 

r 2 = e^-^^cot/3 
sin a u 



3 For alternative neutrino mass models with 126 representations see [17]. 



= h + f (6) 
= hr x + fr 2 
= hri — 3r2/ 
= h-3f 
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The mass sumrules in Eq. (7) provide the first important ingredient in discussing the 
neutrino sector. In the case without any phases in the Yukawa sector, they determine 
completely the input parameters of the model. 

To see this let us note that Eq. (7) leads to the following sumrule involving the charged 
lepton, up and down quark masses: 



where k and r are complex numbers which are functions of the symmetry breaking param- 
eters of the model; the mass matrices M Uy d,i are general symmetric complex matrices. In 
the Eq. (8), tilde denotes the fact that we have made the mass matrices dimensionless by 
dividing them by the heaviest mass of the species i.e. up quark mass matrix by m t , down 
quark mass matrix by m b etc. 

We now proceed to do the phase counting in the model. First we absorb the phase of r 
by redefining k, so that r becomes real. We then choose a basis so that Md is diagonal and 
real. The (u, c, t) basis is appropriately defined so that the weak current is diagonal in this 
basis and M u is still a general complex symmetric matrix wherein the CKM mixing matrix 
is buried. We can now write M u = V T M^ mg V where V is a general 3x3 unitary matrix, 
which has three real rotation angles and six phases. Three of these phases can be put into 
the diagonal elements of the down quark mass matrix and two can be put into the three 
diagonal elements of M u and one remaining phase is in the CKM matrix. The V matrix can 
now be parameterized as V = PdUcKuPu, where P u> d are diagonal unitary matrices. 
Input parameters in the model: 

• We have two parameters (k, r) and six phases constrained by the fact that they must 
reproduce the correct charged lepton masses and lead to large neutrino mixings via 
b — t mass convergence. 

• The other degrees of freedom arise from the fact that the quark masses and mixings 
extrapolated to the GUT scale have uncertainties in them. 

• We have the freedom to change the sign of the quark and lepton masses, which amount 
to to a redefinition of the fermion fields by a 75 transformation. 

We use the above input parameters to get the correct charged lepton masses and subse- 
quently via the type II seesaw (explained below) to predict the neutrino mixings. As it is 
apparent, this is a highly nontrivial task and restricts the parameter space of the theory i.e. 
values of quark and charged lepton masses at the GUT scale, very strongly. 

III. TYPE II SEESAW FORMULA AND MAXIMAL NEUTRINO MIXINGS 

In this section, we explain the type II seesaw formula that we use in discussing neutrino 
mixings. The familiar seesaw formula (type I seesaw) for small neutrino masses [4] is given 



kMi = rM d + M, 



u 



(8) 



by 



M u = -M v oM^Mt 



T 



(9) 



,D 



8 



where M^ R = fvBL, Vbl being the scale of local B-L symmetry breaking. On the other 
hand, it was pointed out in 1980 that in theories with asymptotic parity conservation, the 
seesaw formula has an additional contribution [18] i.e. 

M v ^ fv L - M vD {Jv BL y x M T v ». (10) 



V 2 



where vl = Note that the / matrix is common to both terms. It is possible to find 

different regions of parameter space where the first or the second term may dominate. For 
example, it was shown in Ref. [23] that when parity symmetry is broken at a much higher 
scale than the SU(2) R symmetry (e.g. by breaking SO (10) via 210 Higgs field), one recovers 
the type I seesaw formula. 

In this section, we would like to make two points: first, the origin of the first term in the 
SO (10) theory under consideration and secondly, conditions under which the first term may 
dominate the neutrino mass. 

In our model, the mass of the right handed neutrino comes from a renormalizable term 
of the form v c v c A\ where A° R is the neutral member of an SU(2)r triplet; parity symmetry 
of the theory then implies that there must be a coupling in the theory of the theory of the 
left handed neutrino of the form vvA° L) where A° L is the neutral member of the left handed 
triplet. If the A° L has a vev, then we get type II seesaw formula. 

To show that in our model, A° L has a vev, let us look at the gauge invariant Higgs field 
terms in the superpotential. First we note the decomposition of the 126 under the group 
SU(2) L x SU(2) R x SU{A) C : 



126 = (1, 1, 6) © (2, 2, 15) © (3, 1, 10) © (1, 3, 10) (11) 

The SU{2)l triplet that contributes to the type II seesaw formula is contained in the mul- 
tiplet A L = (3,1,10) and it couples to the left handed multiplet ip = (2,1,4) of the 16 
dimensional SO(10) spinor that contains the matter fermions i.e. i/jliPlAl- On the other 
hand the mass of the RH neutrinos comes from the coupling of A^ = (1,3, 10) submultiplet 
of 126 to the right handed lepton doublets. 

The vev of the neutral member of Ar breaks the B-L symmetry and gives mass to 
the RH neutrinos. This generates the second term in the type II seesaw formula. To 
see how the A° L vev arises, note that the general superpotential of the model contains 
terms of type Ail26 2 • 54 and A 2 10 ■ 10 ■ 54. In the Higgs potential, this generates a 
term (from I-F54I 2 ) of the form 10 • 10 • 126 ■ 126. In this expression, there is a term of the 
form 0(2, 2, 1) 2 A^(3, 1, 10)A^(1, 3, 10) with a coefficient A1A2. Furthermore, in the Higgs 
potential, there is a mass term for A^(3, 1, 10) of the form n\ + X^vfj, where vjj is the GUT 
scale. On minimizing the potential, these two terms lead to a vev for the SU{2)l triplet 

«■ -< *l >- 

It is now clear that if we choose A 3 such that n\ + X 3 v fj <C v BL , then the entries in the 
second matrix in the type II seesaw formula can much smaller than the first term. When 
this happens, then Eq. (6) can be used to derive the sumrule 

Ml = a(M e - M d ) (12) 

This equation is key to our discussion of the neutrino masses and mixings. 
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A. Maximal neutrino mixings from type II seesaw 



Using Eq. (12) in second and third generation sector, one can understand the results of 
[19] in a heuristic manner as follows. The known hierarchical structure of quark and lepton 
masses as well as the known small mixings for quarks suggest that the matrices for the 
second and third generation have the following pattern 



A 2 A 2 

2 




(13) 



where A ~ 0.22 (the Cabibbo angle) as is required by low energy observations. It is well 
known that in supersymmetric theories, when low energy quark and lepton masses are 
extrapolated to the GUT scale, one gets approximately that m b ~ m T . One then sees from 
the above sumrule for neutrino masses Eq. (12) that all entries for the neutrino mass matrix 
are of same order A 2 leading very naturally to the atmospheric mixing angle to be large. 
Thus one has a natural understanding of the large atmospheric neutrino mixing angle. No 
extra symmetries are assumed for this purpose. 

For this model to be a viable one for three generations, one must show that the same 
b — T mass convergence at GUT scale also explains the large solar angle 9 12 and a small 9 13 . 
This has been demonstrated in a recent paper [20]. 

To see how this comes about, let us ignore the CP violating phases and recall Eq. (8). 
Note that in the basis where the down quark mass matrix is diagonal, all the quark mixing 
effects are then in the up quark mass matrix i.e. M u = U^ km M^Uckm- Using the 
Wolfenstein parametrization for quark mixings, we can conclude that that we have 

/A 4 A 5 A 3 \ 



M d « m b 



(14) 



A 5 A 2 A 2 
\X 3 A 2 1 J 

and Mi and M d have roughly similar pattern due to the sum rule 8. In the above equation, 
the matrix elements are supposed to give only the approximate order of magnitude. As we 
extrapolate the quark masses to the GUT scale, due to the fact that — m T m r A 2 for 
some value of tan f3, the neutrino mass matrix M v = ciMd — M() takes roughly the form 
[24] 



A 4 


A 5 


A 3 


A 5 


A 2 


A 2 


A 3 


A 2 


A 2 



M v = c(M d - M e ) « m A 5 A 2 A 2 (15) 

\a 3 a 2 a 2 / 

It is then easy to see from this mass matrix that both the 9 12 (solar angle) and 9 2 s (the 
atmospheric angle) are large. Furthermore, if only terms of A 2 are kept, then, m 3 rs cA 2 and 
in the limit of maximal atmospheric mixing and both m 2 and m\ vanish. As soon as terms 
of order A 3 , ra^ pick up mass and then one has m 2 /m 3 ~ A. This then naturally explains 
the milder hierarchy among neutrinos compared to that among quarks. 

The detailed predictions of the model such as the magnitudes of these angles and neutrino 
masses m 12 depend on the details of the quark masses at the GUT scale and we discuss it 
in the following sections for different cases. 
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IV. PREDICTIONS FOR NEUTRINO MIXINGS WITHOUT ANY CP PHASES 

IN THE MASS SUMRULES 



Let us first consider the case where the Yukawa couplings in the superpotential and the 
vevs of doublet Higgs are all real. This case was considered in our earlier paper [20]. In this 
case all CP phases needed for understanding the observed kaon and B CP violation arise 
from the supersymmetry breaking sector. We can start by solving for the parameters k and 
r in Eq. (8) and find the range of quark masses for which the charged lepton masses come 
out right. We then use the values of k and r as well as the quark masses to get the neutrino 
masses and mixings using Eq. (12). 

While all our predictions are done via detailed numerical analysis using Mathematica, in 
this section we provide a qualitative discussion of the nature of the constraints on the model 
parameters. The qualitative discussion brings out several things clearly: 

(i) while the masses of tau lepton and muon fix the values of k and r, getting electron 
mass is nontrivial and requires fine tuning for quark masses and also the V uh parameter 
within the range allowed by present data. 

(ii) Secondly, we derive an approximate form for the neutrino mass matrix and and show 
how the model generally tends to predict values of U e j, close to its present upper limit. 

To find k and r numerically, we need to specify the values of the quark and lepton masses 
as well as the CKM mixings extrapolated to the GUT scale. These have been discussed 
extensively in literature. We use the values from the paper of Das and Parida [25] and are 
given in Table. 



input observable 


tan (3= 10 


tan f3 = 55 


m u (MeV) 


79+ u - 13 

u - ' z -0.14 


u - ' z -0.14 


m c (MeV) 


210.32™ 


210.50™ 


m t {GeV) 


82.43lSre 


95.14™ 


m d (MeV) 


i - ou -0.23 




m s (MeV) 


29.94^54 


29.81_4' 4 g 


m b (GeV) 




1 41+U.48 


m e (MeV) 


0.3585 


0.3565 


m M (MeV) 


75 671 ^+°- U57 « 


75.29381& 


m T (GeV) 


1 OQ99+0.U013 
1 - zyzz -0.0012 


1 R9Q9+U-U443 
1.UZ»Z_q 0294 



Taking the 2-3 submatrix of the Eq. (8) and remembering that all mixing angles are small, 
we get from this equation the approximate equations: 

k~r + l (16) 

m u m s m c 
k — - ~ r 1 

m T nib m t 

Using these equations and rough numbers for m^/rrir ~ 0.059 and m s /mb ~ —0.026, we get, 
k ~ 0.29 and r ~ —0.71. The results of our detailed numerical solutions for this case give 
-0.78 < r < -0.74 and 0.23 < k < 0.26. 
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Now we illustrate with a particular choice of fermion masses how values of k and r are 
determined and then show how the small electron mass comes out in the model. For this 
purpose, let us express the Mg in terms of the small parameter A ~ 0.22 (the Cabibbo angle) 
for a particular choice of m b (Mu) ~ 1.2 and m T ~ 1.29 as an example: 

^1.6A 4 ±11A 6 

Mi = | 6A 5 T2.1X 2 3.3A 2 | (17) 

4.5A 3 

where we have kept the ± sign in M^n to reflect the sign freedom in the quark masses. 

An approximate expression for the electron mass fh e can be obtained from the above 
expression to be: 

S/T M? 2 13 /10 x 

m e ~ M ttU - — — (18) 

~ Tl£X 4 sgn{m d ) ± sgn(m c )l0\ 6 - 20A 6 

Since rh e ~ A 5 , we need a cancellation between the different terms in the above equation. 
Note that if < and m c < 0, then we get cancellation between the three terms in the 
above equation and we can get the right value for the electron mass. It must however be 
stressed that m e is reproduced only for particular choice of the bottom and tau masses at 
the GUT scale and of course these are in the allowed domain but nonetheless they reflect 
the constraint on the mass parameters for the model to be acceptable. It is impressive that 
it works. 

We caution that the above discussion is meant to give a flavor of the constraints on 
the model. In the detailed numerical calculations, the range of parameters where correct 
electron mass results is larger than what would be implied by the above discussion. 

To study the neutrino masses and mixings, let us write down the neutrino mass matrix 
in this model: 



M v 



/ym d -z(m c X 2 + A 2 X (i \A\ 2 ••• •• \ 

zm c \ + zA 2 \ 5 \A\ ym s — zm c — zA 2 \ 4 '■■ 

\ zAm c \ 3 + zA\ 3 \A\ -zA\ 2 (l + fn c ) y - zrh c A 2 A 4 - z J 



(19) 



where Ay = {m b /m T — r/k) and 4z = k 1 and A = (1 — p — ir/). Noting that m d ~ A 4 , 
m c ~ A 4 , the neutrino mass matrix to leading order takes the form 

/ zA\ 3 \A\\ 

M v = ym s zA\ 2 (20) 

\^A 3 |A| zA\ 2 y-z J 

From this we see that maximal neutrino mixing as well as the correct mass hierarchy comes 
out. Obviously to get the mixing angles in the desired range, detailed analysis is needed 
and we have carried it out. 

Secondly, one can give an "analytic" argument that U e 3 will be close to its present upper 
limit. Again Eq. (8) comes in handy. Roughly U e3 ~ M ^ ;{ h ]i ^ ^(mJ^m^) - A - 



12 



A. Detailed Numerical Analysis 



We undertake extensive scanning of the parameter space of the model defined by the 
uncertainties in the values of the fermion masses at the GUT scale and the values of k and 
r in the neighborhood of the values given above. We have used the values of the standard 
model fermion masses from Table I and the following values for the mixing angles. 



UcKM 



f 0.974836 0.222899 -0.00319129 \ 

-0.222638 0.974217 0.0365224 
V 0.0112498 -0.0348928 0.999328 j 



(21) 



The strategy in our numerical calculations is the following: we focus on the Eq. (8) and 
using the inputs in the right hand side, we look for the eigenvalues of the charged lepton 
mass matrix in the left hand side to match the observed lepton masses. Since there are only 
two unknowns, the second and third generation masses largely fix the parameters k, r as we 
just described. To match the electron mass is nontrivial. The parameters are essentially the 
signs of the fermion masses and the present uncertainties in the values of the SM fermion 
masses. After we fix the electron mass it narrows down our parameter range somewhat. We 
then look for neutrino masses and mixing angles using Eq. (12) in the remaining (small) 
parameter range. Note that due to overall scale freedom of the type II seesaw scale, we 
cannot predict the Am 2 A . We also do a direct numerical solution. Both the results are in 
agreement. 

The solutions we present here correspond to m e ^ T ^ t > and m c ^ s < up to an overall 
sign. 

Our results are displayed in Fig. 1-3 for the case of the supersymmetry parameter 
tan j3 = 10. In these figures, we have restricted ourselves to the range of quark masses for 
which the atmospheric mixing angle sin 2 29a > 0.8. (For presently preferred range of values 
of sin 2 29a from experiments, see [26]). We then present the predictions for sin 2 2# , Am 2 
and C/ e 3 for the allowed range sin 2 29 A in Fig.l, 2 and 3 respectively. The spread in the 
predictions come from uncertainties in the s, c and the 6-quark masses. Note two important 
predictions: (i) sin 2 29 Q > 0.91 and U e ^ ~ ±0.16. The present allowed range for the solar 
mixing angle is 0.7 < sin 2 29 Q < 0.99 at 3er level [26,27]. The solutions for the neutrino 
mixing angles are sensitive to the b quark mass. 

It is important to note that this model predicts the U e s value very close to the present 
experimentally allowed upper limit and can therefore be tested in the planned long base line 
experiments which are expected to probe U e3 down to the level of ~ 0.05 [28,29]. 



V. EFFECT OF CP PHASES ON NEUTRINO MIXINGS 

In this section, we keep the CP phases as described before and look for solutions to the 
charged lepton equation, Eq. (8) and then use the allowed parameter range to look at the 
predictions of the neutrino masses and mixings. We have seven CP phases including the 
CKM phase. One might think that since there are more parameters in this case, getting a 
solution will be trivial. We find that actually, it is not so. Let us demonstrate this in an 
analytic way very crudely and we follow it up with detailed numerical scan to get the mixing 
angles. 
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□ . 82 0.84 . 86 0.88 0.9 

3in Z (26 A ) 

FIG. 1. The figure shows the range of predictions for sin 2 29q and sin 2 29 a for the range of 
quark masses in table I that fit the charged lepton spectrum and where all CP phases are set to 
zero. We required that the ratio Am^/Am^ < 0.05. Note that sin 2 29 & > 0.9 and sin 2 29a < 0.9. 



To proceed with this analysis, first note that the phases are distributed in Eq. 
(8) as follows. First the down quark mass eigenvalues are chosen to be complex i.e. 
(m^e*^ 1 , m s e J/32 , e l133 ) and similarly for the up quark sector, there are three phases in the 
masses (denoted by ctj). Using 8, we see that 



kMf 



( rm d e iPl + m c X 2 e ia2 m c Xe ia2 + A 2 X 5 e ias (A) Am c X 3 e ia2 + AX 3 e ia3 (A) \ 
+A 2 X 6 e ia *(A) 2 

rh c \e ia2 + A 2 X 5 e ioi3 (A) rm s e if32 + m c e ia2 + A 2 X A e ia * m c AX 2 e ia2 + AX 2 e ia3 

V Am c X 3 e ia2 + AX 3 e ia3 (A) rh c AX 2 e ia2 + AX 2 e ia3 re iA + m c A 2 X 4 e ia2 + e ias J 



where we have used the Wolfenstein parametrization for the CKM matrix. Essentially the 
eigenvalues of the matrix on the right hand side must match the charged lepton masses. 
Note that as in the case without the CP phase, this will require us to work only in a very 
limited range of the quark masses. Since now we have to align the phases with the known 
standard model CKM phase, the constraints are even tighter. The muon and tau masses 
come out quite easily as in the case without phases discussed in the previous section. 

As in the CP conserving case, the mass of the electron requires fine tuning as can be 
seen by noting that the 11 entry of Mi is roughly of order A 4 which is about 3 to 4 times 
the value required to give the correct electron mass. One must therefore have cancellation 



(22) 
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0.S2 0.84 0.86 tl.88 0.9 

3in Z (26 A ) 

FIG. 2. The figure shows the predictions for sin 2 29 a and Ara^/Ara^ for the range of quark 
masses and mixings that fit charged lepton masses and where all CP phases in the fermion masses 
are set to zero. 



to get the correct m e . Before discussing the details of this cancellation, let us first look at 
the neutrino mass matrix, which looks as follows: 



I zAX 3 {l - p - ir])e io 

0.7y\ 2 e^ zA\ 2 e ias j (23) 

v zA\ 3 (l - p - ir])e ia3 zA\ 2 e iai ye 1 ^ - ze ia * 



First thing to note is that we must have 0:3 = ^3 for b — r unification to lead to maximal 
mixing. 

To see very qualitatively what value of CKM phase is required for the model to work, 
we need to analyze the charged lepton masses. As in the CP conserving case, the muon and 
the tau mass come out for the choice of k « 0.25 and r ps —.75. Let us now discuss the 
electron mass. For this we first note that we use the Wolfenstein parametrization. The two 
parameters p and rj responsible for the CKM phase are given by the latest data to be [21] 
.12 < p < .35 and 0.28 < i] < .41. Since the charged lepton mass matrix has hierarchical 
structure, we can deduce an approximate formula for m e from the sum rule 8 and using the 
above values of k and r: 



~ -3m d e l/31 + m c \ 2 e ia ' 2 + A 2 \ 6 e ia3 (A) 2 - (4AA 6 (1 - p - ir]) 2 e ia3 ) - 0(A 7 ) (24) 



Numerically, to get the correct value for m e ~ 0.00028 ~ 0.6A 5 , we must cancel most of the 
rh d contribution to it i.e. Srhd — 1.3A 4 , from the other terms: these terms get maximized 
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-0. 15 



0.86 

Sin" (2e A ) 

FIG. 3. The figure shows the predictions of the model for sin 2 29a and C/ e 3 for the allowed 
range of parameters in the model as in Fig. 1 and 2. Note that U e s is very close to the upper limit 
allowed by the existing reactor experiments. 



for p < in which case we get |1 — p — irj\ ~ 1.41 so that the last term contributes 
32A 6 ~ 1.5A 4 . It also determines the arbitrary phases that accompany the quark masses 
as follows: j3\ = «2 = «3 + 25 where tan 5 ~ -^— p and < 0. The terms in the electron 
mass formula then cancel leading to rh e ~ .9A 5 which is in the rough "neighborhood" of 
the observed value (rh e ). This value of p gives a CKM phase which is in the third quadrant 
implying that in order to understand the observed CP violation i.e. sin 2(3 in B-system, one 
will need to invoke CP violation from the supersymmetry breaking scalar masses. 

Again we caution the reader that this is a crude estimate. In detailed numerical analysis, 
we can get solutions for small positive p and one next has to see how these values fit the 
neutrino mixings. Thus combined fit to both electron mass and neutrino mixings works only 
for negative p value which is different from the current standard model CKM fit. 

The detailed numerical predictions for this model are given in Fig. 4,5 and 6, where we 
have scanned over all the phases looking for the correct charged lepton masses and acceptable 
neutrino parameters. We have also extended the domain of Amg/Am^ to 0.08. 

In concluding this section, we note that since the main constraints on the model came 
from trying to fit the electron mass, it is worth pointing out that one way to avoid these 
constraints is to include the effects of the higher dimensional operators (HDO) which con- 
tribute to the electron mass and is too small to be of relevance in the discussion of other 
masses. Since usually the contribution of HDO terms in this theory would contribute terms 
of order (-^-J — 10 _2n , where 3 + n denotes the dimensionality of the HDO. The leading 
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FIG. 4. The figure shows the predictions for sin 2 29 a and sin 2 29® for the range of quark masses 
and mixings that fit charged lepton masses in the presence of all CP phases in the fermion masses. 
The four different panels give the predictions for different values of the CKM phase (tan 5c km = «)■ 
Note that all these values are outside the one sigma region of the present standard model fit to all 
CP violating data. Note that the case Sckm = 180 includes the CP conserving case discussed in 
sec. IV. 



order HDO must therefore correspond to n — 2. This can be guaranteed for instance by 
demanding that superpotential is R-odd and all the Higgs fields are R-odd, whereas the 
matter fields are R-even. We have repeated our computations with the HDO terms men- 
tioned and we find that one can obtain the neutrino mixing parameters in the right range 
for the standard model value of the CKM phase. 



A. CP violation in leptonic mixings 

The neutrino mass matrix discussed above can be approximately diagonalized to obtain 
the CP phases in the PMNS matrix. For this purpose, we first diagonalize the charged 
lepton mass matrix. The mixing matrix has a form similar to the CKM matrix. Therefore, 
in the basis where charged leptons are mass eigenstates, the neutrino mass matrix remains 
approximately of the same form as in Eq. (23) i.e. 



M u = m zAX 2 



( ae i(S+a 3 ) 

(l + e)e iA e ias j (25) 
\ae e zAx2 
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FIG. 5. The figure shows the predictions for sin 2 26 a and Am^/Ara^ for the range of quark 
masses and mixings that fit charged lepton masses and where all CP phases in the fermion masses 
are kept subject to the condition that b-tau mass convergence is responsible for large neutrino 
mixings, for four different values of the Sckm- 



where a 



l±p) 2 + V 2 zAX; 



e ~ — 1 and tan 5 ~ . In the expression for a 
and tan S, the plus sign corresponds to the case without any contribution from the higher 
dimension terms to the electron mass and the minus corresponds to the case with higher 
dimensional terms. In order to get maximal mixing out of b-tau mass convergence, we must 
have a 3 = @ 3 and (ye' l/33 — ze' ias ) ~ zAX 2 . Furthermore, to get the solar mass difference right, 
we must have /3 2 = ct 3 . We can therefore take the common phase e %az out of the matrix 
which leaves us with a simple matrix of the form 



M v = 



I 


\ ae 



iS 





1 



ae 
1 
1 



iS 



J 



diagonalize. The matrix that diagonalizes it is of the form 



U 



V 



V 



c — s 

( S (l-t)W2|) (-c(l-f)W2|) 
(s(l+j)- c> /2f) (c(l+f)W2|) 



V2 



V2 



(26) 



IV2\ 

j_ / 

V2 ) 



(27) 



where V = diag(e , 1,1). Note that since the PMNS mixing matrix is real to order A, 
direct CP violating effects among neutrinos are suppressed and unobservable in near future. 
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FIG. 6. The figure shows the predictions for sin 2 29a and U e 3 for the range of quark masses 
and mixings that fit charged lepton masses and where all CP phases in the fermion masses are 
included, for four different values of 5ckm- 



This is a definite prediction of the minimal SO (10) model. In a sense this is not surprising, 
since the lepton sector and the quark sector are linked by quark-lepton unification and we 
know that in the quark sector the CP violating phase is only present in the A 3 order. 

Since the model predicts a large U e ^, if there is any CP violating phase in the mixing 
matrix, it has a better chance of being observable; on the other hand this model predicts 
no CP violation, so that the model has a better chance of being falsifiable and strengthens 
the case for searching for CP violation in neutrino oscillations. It is important to emphasize 
that the model however enough CP violation to in the right handed neutrino couplings so 
that one should expect enough baryogenesis. 



VI. SUMMARY AND CONCLUSIONS 

In summary, we have shown that the minimal SO(10) model with fermion masses re- 
ceiving contributions from only one 10 and one 126 Higgs multiplets is fully predictive in 
the neutrino sector. It predicts all three neutrino mixing angles in agreement with present 
data but with a value of U e z which is very close to the present upper limits from the reactor 
experiments. This high value of U e3 provides a test of the model. We also find that the in- 
troduction of CP phases in the Yukawa couplings still keeps the model predictive. The CKM 
phase in this case is outside the one a region of the present central value in the standard 
model suggesting that there are new CP violating contributions from the SUSY breaking 
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sector. We find it interesting that to the leading order in the Cabibbo angle, the leptonic 
CP violation vanishes in spite of the fact that the quark sector has CP violation. 

There are no additional global symmetries assumed in the analysis. The neutrino data 
once refined would therefore provide a crucial test of minimal SUSY SO (10) in the same 
way as proton decay was considered a crucial test of minimal SU(5). 

This work is supported by the National Science Foundation Grant No. PHY-0099544. 
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